Evolution of a Bose-Einstein condensate of neutral atoms 
— A field theoretical approach 
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The particle distribution in a Bose condensate under the 
trapping potential and its time evolution after switching off 
the trapping potential suddenly are calculated. We investi- 
gate the problem from the viewpoint of quantum field the- 
ory,@using a model of self-interacting neutral boson field. 
Within the approximation of retaining the most dominant 
term in the Hamiltonian after applying the Bogoliubov re- 
placement, we can calculate analytically the particle distribu- 
tion as a function of space and time coordinates. 



The recent observations of Bose-Einstein condensation 
(BEC) in magnetically trapped atomic gases has 
stimulated the theoretical study of this phenomena. In 
this work we investigate the particle distribution in a con- 
densate of neutral boson field under the trapping poten- 
tial and its time evolution immediately after the trapping 
potential is switched off, taking the interaction between 
the particles into account. 

We consider a collection of interacting Bose ^particles 
trapped in a harmonic oscillator potential | 2i=i ^f^f » 




g / dxdx' cj)\x)(j)\x')V{x - x')<t>(x)(t>{x'), (1) 



whose Hamiltonian is given by 



where we set particle mass= fi = 1. 

In a classical paper on BEC || , Bogoliubov treated a 
similar Hamiltonian but without the trapping potential. 
The key idea in his theory is that the creation and anni- 
hilation operators for the zero momentum particle should 
be replaced by a c-number y/No, when BEC occurs. The 
quantity Ao is the average number of particles occupy- 
ing the zero momentum state and BEC means that this 
number becomes very large and macroscopic. Roughly 
speaking, this replacement is justified as follows: For 
BEC state, A particles are in the zero momentum state. 
Thus the creation and annihilation operators for zero mo- 
mentum state give 



cl\N > = v^VoTTl No + 1 > 
c |A > = v/AoIAq - 1 > . 



(2) 



Since the number Nq is very large, Ao + 1 ~ Ao. There- 
fore both the operations of cj and Co on system states are 
similar, giving the result y/ Ao as an eigenvalue. Thus we 
can treat both the operators cj and cq as c-number \ZNq. 

To extend this idea to a trapped Bose gas, first let 
us introduce annihilation operators of particles a n in the 
harmonic oscillator modes by expanding the field oper- 
ator 4>(x) in terms of a complete set of the harmonic 
oscillator wave functions u n (x): 

(j>{x) = y^ j a n u n {x), (3) 

n 

where u n (x) satisfies 
1 1 3 

{~2 A+ ^E^iKW =e n u n (x), (4) 
i=i 

and n = (ni,ri2,n>3) with non-negative integers rij. 

In a condensed state, most particles are in the ground 
state of the harmonic oscillator, that is, the state for n — 
(0,0,0). Therefore, the occupation number Ao of this 
state is supposed to be very large. In typical experiments, 
the total particle number is of order of 10 3 ~ 10 5 . Then 
we can apply the Bogoliubov replacement to the creation 
and annihilation operators, a and oq : 



Then the field operator Mx) reads 



x) = y/Nou (x) +y^/a n u n (x), 



(5) 



(6) 



where the prime in ^ n ' means to omit the term n — 
(0,0,0). 

Substituting the above expression of 4>(x) into Hb in 
(1), we have 

H B =N e + ^2 'snaicin 



1 



+ \ aN Q J dxdx ' u o( x ) v ( x - x ') u o( x ') 
+gN 3 Q /2 J2'(ai + a n ) 



x J dxdx' u n (x)uo(x)V (x — x')v,q(x') 

+O(gN ) (7) 

when V(x) = V(—x) is assumed. Since No is large, we 
shall retain only terms of order Nq and N^ 2 in the in- 
teraction part. We remark the presence of the terms of 

3 /2 

order N ' .In the original Bogoliubov theory, the terms 

3/2 

of order N do not exist, which is due to the momen- 
tum conservation. In the present case, however, we do 
not have such a conservation law, and the dominant q- 
number terms are linear in the creation and annihilation 
operators of excited states. 

By defining operators cv n and a n by 



No^n 



(8) 



which satisfy the same commutation relations as and 
a n , we can diagonalize the approximated Hamiltonian 
thus obtained: 



H b = 'e n oi\oL n + 7V eo(l + ^7o 



with the c-number coefficient 

7 



In terms of these a's, the particle number density op- 
erator h(x) is written as 

n(x)=(fy (x)(/)(x) 

=^ 'al l a m u n (x)u rn .(x) 

+ ^/No~ 'i a l + a n )u n (x) 

n 

X |u (x) - y^TrnMrnQE) j 
+N (u (x) - J2 'fnU„(x^J (11) 

and its thermal average, 

n(x) = Z~ lr Fr [e~ pHB n{x)\ (12) 

with Z — Tr [e - ^- 8 ] , can easily be calculated by making 
use of (§) to find 

n(x) = N Q ^uo(x) -^2'~/ n u n (x)^j 

+ E'^tt ( 13 ) 



In our approximation of taking only the linear terms in 
a' n and a n , the energy spectrum of single particle states 
are not affected by the existence of interactions between 
particles but the wave function of the ground state does 
change. This can be seen in the first term above. 

The total particle number is also obtained by integrat- 
ing the above n(x) with respect to x: 



N = No(l + J2'ri)+J2' 



1 



(14) 



This equation gives the relation between the temperature 
and the condensate fraction /o = No/N. It shows that 
the interaction between particles lowers the fraction fo- 
Now we shall investigate how n(x) evolves, if the trap- 
ping potential is removed. Let us switch off the trapping 
potential at t = 0, thereafter the system develops under 
the Hamiltonian, 



H=jdx^(x)S.-^cf>(x) 



+ ^g J dxdx' ${x)<j)\x')V{x - x')(t>{x)(t>{x'). (15) 

For a while, however, the system must be still in BEC. 
Therefore we may be able to approximate the self- 
interaction term (the second term) in ( |l5| ) by one in (Q) 
(to neglect the terms of order gNo). Hence we have 



■„ = ^-J dxdx'u n (x)u (x)V(x - x')ul{x'). (10) H - j 



lh, I dx^(x){-^A}<j)(x) 

No^2'(al + a n )e nln + N Q eo(l + \lo)- (16) 



Here Nq and hence 7„ must become time dependent. 
In principle, this dependence will be determined self- 
consistently by the condition that the total particle num- 
ber is conserved. In the following, however, our consider- 
ations are limited to the case of time- independent No and 
7„, because we are interested in the short time evolution 
after switching off the trapping potential. 

The quantum number n does not diagonalize the ki- 
netic term in (|l5|). It is convenient to introduce the 
creation and annihilation operators in momentum eigen- 
states cj, and Ck in a usual way: 



0(a) 



Inversely we get 



Ck 



1 



(2tt) 3 / 2 



dk Cke 



k x 



(17) 



dxe- ikx 6(x) 



(2tt) 3 /2 

(k), 



(18) 



where we use the Bogoliubov replaced field operator 
The operators a n and Ck are related as, 
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a n = dxu n (x)4>(x) = dkc k u* n (k) 



where u means Fourier transform of u denned by 



1 



(2tt) 3 / 2 



dx f(x)e~ 



ik x 



(19) 



(20) 



Putting @ and (Q) into (|T|), we have an expression 
of H written by c£, and c k as 

+iV eo(l + ^7o) (21) 

with £fe = fc 2 /2. 

The Heisenberg operator for t > 0, 

c fe (i)=e^c fe e- itif (22) 

satisfies the equation of motion 

ic k {t)=e UH [c k ,H]e- UH 

(k). (23) 

n 

We can solve this equation with the condition c k (0) = c k 
to get 



c/t(i) = c fc e- fe ** + 



1 



^J^'e^nM*)- (24) 



Substituting ( |18| ) into the first term above, we have 

n 

e —ie k t _ ^ 



/ jV Q ^ yg n 7nMn(fc). 



(25) 



The field operator in the Heisenberg picture is written as 



(x,*)=e rt -" <f>(x)e 



-itH 



1 



e?fc c k (t)e 



k x 



(2tt) 3 / 2 

=^ 'a n u n (x, t) + ^/J%u (x, t) 

-iy/flo y^'e n 7n / d\u n (x,\), 
Jo 



where we denote 



u n (x,t) 



1 



dfe u n (k)e 



kx 



(2tt) 3 / 2 , 

In terms of a„'s, (Bq) can be rewritten as 



(26) 



(27) 



<j>[x, t)=^2 'a n Un(x, t) + ^/Nou n (x, t) 

n 

-y/No^2'"tnU n {x,t) 

n 

-iy/No5^'e„n n [ d\u n {x,\) (28) 

and the number density operator in the Heisenberg pic- 
ture is given by 

n(x, t)=(jr(x, t)(f>(x, t) 

'ala m u^(x, t)u m (x, t) 

71 

+N \f(x,t)\ 2 , (29) 

where 

f{x,t) = u (x,t) 

- y2'^ n {u n (x,t) +ie n / dXu n (x,X)}. (30) 

The thermal average of n(x, t) is given by 

n(x, t) = Z~ lr Tr [e~ mB h(x, t)] , (31) 

since the density matrix does not change in time in the 
Heisenberg picture. Making use of the eqs.(||) and (ps|), 
finally we find 



n(x, t) = 
N \u {x,t) - 
1 



E 



'j n {u n {x,t) +ie n / dXu n (x,\)}\ 



\u n (x,t)\' 



(32) 



As seen in the definition (|27j), the functions u n (x,t) 
represent the wave functions of free motion started from 
the harmonic oscillator eigenfunctions. Comparing the 
above n(x,t) with n(x) given in ([l3"|), we can see that 
n(x,t) is the result of free evolution of n(x), except for 
the term ie n f*dXu n (x 7 A). This term gives the effect of 
the interaction between particles. If we drop this term, 
the integration of n(x, t) with respect to x becomes the 
same as N in eq.(|14|), because of the orthogonality of 
{u n (x,t)}: 



dxu* n (x,t)u rn {x,t) = 5„, 



(33) 



This means that, in our approximation, the time depen- 
dence of the fraction /o stems from the interaction be- 
tween particles. 

In summary, the two main results of this work are as 
follows: (I) The ground state particle distribution in a 



3 



Bose condensate under the trapping potential, seen from 
the first term in (|l3| ) , is different from the noninteracting 
case, where the distribution is expected to be |ito(a;)| 2 . 
(II) The time evolution of the particle distribution of a 
condensed Bose gas immediately after switching off the 
trapping potential suddenly can be calculated analyti- 
cally as in (p^). In ([l3|) and (^|) the parameter j n , rep- 
resenting the interaction effect, is crucial quantitatively. 
Simple manipulations of (|l0|) lead to 



for all even rij 
= otherwise (34) 

where g = Aira and a is the s-wave scattering length. To 
estimate the numerical order of the interaction effects in 
our calculations, let us see the values of the parameters 
in the experiment [jjj using Rb atoms. There we have 
No = 2 x 10 3 , a = llOao (ao:Bohr radius), loz = %/8wi = 
V&J2 = 2007rHz, then the leading 7„ for n — (2, 0, 0) or 
(0, 2, 0) and the next-leading one for n = (2, 2, 0) become 
-0.23 and 0.06, respectively. 

Our calculations are based on Bogoliubov theory, 
adapted to a system trapped by a harmonic oscillator 
potential. Our approximation scheme is to take only the 

— 1/2 

dominant terms in N expansion, but is still nonper- 
turbative with respect to g (the interaction between par- 
tides). Higher order calculations in N expansion is 
under study. The future task is to combine this field 
theoretical approach with nonequilibrium Thermo Field 
Dynamics 

This work was contributed to the International Work- 
shop on Thermal Field Theories and Their Applications 
in Regensburg, 10.-14. Aug. 1998. 
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